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Abstract. We have determined the currently allowed regions of the parameter spaces of two 
representative models of diffuse neutrino flux from active galactic nuclei (AGN): one by Koers 
& Tinyakov (KT) and another by Becker & Biermann (BB). Our observable has been the 
number of upgoing muon-neutrinos expected in the 86-string IceCube detector, after 5 years 
of exposure, in the range 10 5 < E v /GeV < 10 8 . We have used the latest estimated discovery 
potential of the IceCube-86 array at the 5a level to determine the lower boundary of the 
regions, while for the upper boundary we have used either the AMANDA upper bound on 
the neutrino flux or the more recent preliminary upper bound given by the half-completed 
IceCube-40 array (IC40). We have varied the spectral index of the proposed power-law fluxes, 
a, and two parameters of the BB model: the ratio between the boost factors of neutrinos 
and cosmic rays, T v /Tcr, and the maximum redshift of the sources that contribute to the 
cosmic-ray flux, Zqr*. For the KT model, we have considered two scenarios: one in which 
the number density of AGN does not evolve with redshift and another in which it evolves 
strongly, following the star formation rate. Using the IC40 upper bound, we have found that 
the models are visible in IceCube-86 only inside very thin strips of parameter space and that 
both of them are discarded at the preferred value of a = 2.7 obtained from fits to cosmic-ray 
data. Lower values of a, notably the values 2.0 and 2.3 proposed in the literature, fare better. 
In addition, we have analysed the capacity of IceCube-86 to discriminate between the models 
within the small regions of parameter space where both of them give testable predictions. 
Within these regions, discrimination at the 5a level or more is guaranteed. 
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1 Introduction 

Active Galactic Nuclei (AGN) are the most luminous persistent objects in the Universe, 
emitting radiation along almost the entire electromagnetic spectrum, with typical luminosities 
on the order of 10 erg s _1 (see, e.g., [1, 2]). There is evidence that supports the idea 
that AGN are powered by matter accreting onto a central supermassive black hole, with a 
mass between 10 6 and 10 10 times the solar mass [3, 4]. In some cases an enormous amount 
of energy is released in the form of two highly-collimated relativistic jets that emerge in 
opposite directions, perpendicularly to the accretion disc. Although the composition of these 
jets is unknown, it is widely believed that they contain high-energy charged particles, such as 
electrons, protons, and ionised nuclei, which have been accelerated as a result of the repeated 
crossings of the shock fronts that exist within gas clouds moving at relativistic speeds along 
the jets. Such a process would be able to give protons and nuclei energies of up to ~ 10 20 
eV [5, 6]. 

Recently, the Pierre Auger Observatory (PAO) claimed to have detected 69 cosmic-ray 
events with energies above 55 EeV [7] (see also [8]), providing evidence of the anisotropy in 
the arrival directions of utrahigh-energy cosmic rays (UHECRs). Based on the observation of 
29 of these events having an angular separation of less than 3.1° from the positions of AGN 
in the 12th edition Veron-Cetty & Veron catalogue [9], a possible correlation was found with 
AGN lying relatively close, at distances of 75 Mpc or less. Even though the claim on the 
correlation has lost some ground since the first publication of the Auger results [10, 11], it 
still constitutes a possible hint towards identifying AGN as the sources of the highest-energy 
cosmic rays. It is also believed that AGN could be sites of ultra-high-energy (UHE) neutrino 
production. These would be produced in the interactions of UHE charged particles among 
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themselves and with the ambient photons. Therefore, under the assumption that cosmic-ray 
emission is accompanied by neutrino emission [12-14], Auger's claim can be used to normalise 
the neutrino flux predicted by astrophysical models of AGN. 

In the present work, we have focused on two such models of neutrino production that 
take into account Auger's results: one by H. B. J. Koers & P. Tinyakov [15] and another one 
by J. Becker & P. L. Biermann [16], which we will call hereafter the KT and BB models, 
respectively. They differ greatly in their assumptions and, within some regions of their 
parameter spaces, on their predictions of the neutrino fluxes. We have assessed the possibility 
of observing these two fluxes in the km-scale IceCube neutrino telescope at the South Pole, by 
allowing their respective model parameters to vary within given boundaries, and calculating 
the corresponding number of high-energy muon- neutrinos expected in the detector. In doing 
this, we have taken into account the experimental upper bound on the neutrino flux set 
by the AMANDA-II experiment [17], an upper bound set by IceCube [18] in its 40-string 
configuration, and the signal discovery potential of high-energy astrophysical neutrinos in 
the completed 86-string IceCube array. Furthermore, we have also explored the parameter 
space for regions where the event-number predictions from the two production models can 
be distinguished from each other. 

The remaining of the paper is divided as follows. In Section 2, we describe the salient fea- 
tures of the KT and BB models, and show explicitly how the observations from the PAO enter 
the flux normalisation. Section 3 introduces current and envisioned experimental bounds on 
the high-energy extra-terrestrial neutrino flux. In Section 4 we allow the parameters in the 
KT and BB models to vary within given bounds, and calculate the number of muon-neutrinos 
in IceCube predicted by each, while, in Section 5, we present comparative plots of the two 
models in parameter space. We summarise and conclude in Section 6. 



2 Two models of neutrino production at AGN 

AGN have long been presumed to be sites of high-energy neutrino production. In the scenario 
of neutrino production by meson decay, it is assumed that inside the AGN protons are 
accelerated through first-order Fermi shock acceleration [6, 19] and that pions are produced 
in the processes 

p + 7^A+^|^+^ + , n + 7 ^p + ir-, (2.1) 

with branching ratios Br (A + — > pn ) = 2/3 and Br (A + — > nir + ) = 1/3. The neutral pions 
decay into gamma rays through 7r° — > 77, while the charged pions decay into electron- and 
muon-neutrinos through 

tt + -> v M + p + -> + e + + v e + , 7T~ -)■ + ijT -)■ + e~ + T> e + . (2.2) 

The gamma rays thus created may be obscured and dispersed by the medium, and the 
protons will in addition be deviated by extragalactic magnetic fields on their journey to 
Earth. Neutrinos, on the other hand, escape from the production site virtually unaffected by 
interactions with the medium, so that, if their direction could be reconstructed at detection, 
they could point back to their sources. 

If neutrinos are produced by charged pion decay, then, from eq. (2.2), the ratios of the 
different flavours {y x + V x ) to the total flux are 

< : ^ : 4P Vt = 1 : 2 : . (2.3) 



- 2 - 



Under this assumption, by the time neutrinos reach Earth, standard mass-driven neutrino 
oscillations will have distributed the total flux equally among the three flavours so that, at 
detection, 

K ■ K ■■ <hr = 1 : 1 : 1 ■ ( 2 - 4 ) 

New physics effects, such as neutrino decay [20], decoherence [21], or violation of Lorentz 
invariance or of CPT [22-24], could in principle result in large deviations from these ratios. 
In the present work, we have assumed that the ratios at production and detection are given, 
respectively, by their standard values, eqs. (2.3) and (2.4). 

In what follows, we will present in detail two representative models of UHE neutrino 
production at AGN, one by Koers & Tinyakov (KT) and the other by Becker & Biermann 
(BB), both of which make use of the apparent correlation between the directions of UHECRs 
and the positions of known AGN reported by the PAO in order to extrapolate the diffuse 
neutrino flux. 

2.1 Cosmic ray flux normalisation 

The preferred mechanism for cosmic-ray acceleration at AGN is first-order Fermi acceleration 
[6], which results in a power-law differential diffuse cosmic-ray proton spectrum, 



6f (*0 = ^ = Af r E-<* > (2-5) 



with E the cosmic-ray energy at detection on Earth and Ap lS an energy-independent nor- 
malisation constant. The integral of this expression, 



*™ (E th ) = I -^dE ~ Af« (a p - I)" 1 > (2-6) 



it- 



is the integrated cosmic ray flux above a certain threshold energy E t ^. Using experimental 
data, the integrated flux can also be calculated as 

$f S {E th ) = N cvts (E th )/E , (2.7) 

where N evts (E t h) is the number of observed cosmic rays above a given value of -Eth and H is 
the total detector exposure. 

Combining this expression with eq. (2.6) yields for the normalisation constant, 

Af = (Op - 1) Eft- 1 = ^evtsCgp- 1) ^-1 (28) 

We will see in the following two subsections that the relation between the cosmic-ray normal- 
isation constant, Ap lS , and the neutrino normalisation constant, A^, lS , is model-dependent. 

When calculating the proton spectrum from a single point source, we will need to 
weigh the normalisation constant using the detector effective area A that is accessible to the 
observation, which depends on the declination 5 S of the source, i.e., 

Af = *? K " 1) E2~ l - ^K-l) ^-iJ^gl (2 g) 

where we have implicitly defined the integrated flux from a point source, Zp . 
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We will use the latest results from the PAO on the observation of UHECRs [7] to evaluate 
the diffuse and point-source cosmic-ray fluxes. Using data recorded from 1 January 2004 to 
31 December 2009, amounting to an exposure of S = 20370 km 2 yr sr, the total number 
of UHECRs with zenith angles 9 < 60° and reconstructed energies above E t h = 55 EeV is 
Ntot = 69 events. Of these, the arrival directions of iV corr = 29 events were found to lie at an 
angular distance of less than 3.1° from the position of an AGN within 75 Mpc (z < 0.018) 
in the 12th edition Veron-Cetty & Veron (VCV) catalogue, i.e., they were correlated to an 
identified AGN. In particular, A r c C n A = 2 events were correlated to Centaurus A (Cen A), 
the nearest active galaxy, which, at a distance of about 3.5 Mpc, is one of the most promising 
UHE neutrino sources [25, 26]. 

Note that the original PAO report on UHECR anisotropy [8] made use of 9000 km 2 yr 
sr to report a total of 29 events above a threshold energy of 57 EeV, out of which 20 were 
correlated to AGN in the VCV catalogue, and 2 were correlated to Cen A. The neutrino 
production models that we have considered in our analysis were built using these data. In 
what follows, we have updated them using the latest PAO results. 

2.2 Model by Koers & Tinyakov 

The KT model [15] assumes that Cen A is a typical source of UHECRs and neutrinos, and 
computes the diffuse flux under the assumption that all sources are identical to Cen A by 
integrating over a cosmological distribution of sources, while taking into account energy losses 
during the propagation of the particles. Two limiting cases have been considered regarding 
the source distribution: one in which there is no source evolution with redshift, that is, 
e (z) = 1, and another one, adopted from [27], in which there is a strong source evolution 
that follows the star formation rate, i.e., 

r (1 + zf A , if z < 1.9 
e 0) oc \ (1 + 1.9) 3 ' 4 , if 1.9 < z < 3 . (2.10) 
I (z - 3)- a33 , if z > 3 

The integrated UHECR diffuse flux and the integrated flux from Cen A above 55 EeV 
can be calculated, respectively, using eqs. (2.8) and (2.9): 

$f S (E th ) = ^ot-jVcenA = j ^ ^-20 QlR -2 g -l ^-1 (2 n) 

$Cen A {Eth) = NcenAfA (6.) d£l = ^ % ^ ^ ^ 

a A(d s ) 

where 5 S = —43° is the declination of Cen A. The relative exposure at this declination is 
A (5 S ) I f A (5 S ) dQ = 0.15 sr -1 [15, 28]. In eq. (2.11), the number of cosmic rays from Cen A 
is subtracted from the total since the flux is not subject to the energy losses that the diffuse 
flux is, on account of its being the closest AGN. 

The diffuse neutrino flux is normalised using the integrated UHECR flux $f s (E th ) 
above the threshold E^, 

W^W) ~ H iEth) ^enA ( £ th) - (E th ) . (2.13) 

The proportionality constant, H (E^), is called the "neutrino boost factor" and contains the 
information on neutrino mean free path lengths and source evolution. To calculate it, proton 
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energy losses are taken into account in the continuous-loss approximation, considering losses 
by the adiabatic expansion of the Universe and from interactions with the CMB photons 
resulting in pion photoproduction and electron-positron pair production; see Appendix A in 
Ref. [15] for details. The variation of H with a p is shown in the same reference. Note that 
the change in the reconstructed threshold energy from 57 EeV in the original PAO analysis 
[8] to 55 EeV in the updated analysis [7] has reduced H in about 10%. This decrease is 
compensated by a higher value of the ratio &p lS /<&p Cn A , which has moved from 1.8 using the 
original PAO data to 5 using the latest data. As a result, the KT diffuse neutrino flux has 
only changed marginally between the old and new PAO data set. To obtain the diffuse flux, 
the source distribution is integrated up to z = 5. This relation between the diffuse neutrino 
flux and the flux from Cen A is the main result of the KT model. 

In their paper [15], Koers & Tinyakov used a model by Cuoco & Hannestad [29] to 
describe the neutrino emission from Cen A, (f)„ cn A , itself based on a model by Mannheim, 
Protheroe &; Rachen [30]. In this model, it is assumed that high-energy protons, acceler- 
ated by some mechanism (e.g., shock acceleration) are confined within a region close to the 
source. Because of energy losses in their photopion interactions with the ambient photon 
field, which is assumed to have an energy spectrum n{E 1 ) oc E~ 2 , their lifetime is much 
shorter than their diffusive escape time and they decay into neutrons and neutrinos, both 
of which escape the source. Thereafter, the neutrons decay into UHECR protons; however, 
because of their interaction with the photon field before decaying, the neutrons produce a 
softer proton spectrum than the seed proton spectrum. Furthermore, the model predicts two 
spectral breaks in the CR spectrum, at energies at which the optical depths for proton and 
neutrino photopion production become unity. These two breaks are close in energy, though, 
so that, to simplify the model, only one spectral break is considered, at energy E^ r . Below 
Eb r , the UHECR proton and neutrino spectra are harder than the seed proton spectrum by 
one power of the energy, while above E^, the UHECR proton spectrum is softer than the 
seed spectrum by one power of the energy and the neutrino spectrum is harder by one power 
of the energy. Hence, at high energies, the model predicts a neutrino spectrum that is harder 
by one power of the energy than the UHECR proton spectrum. 

Following [15, 29, 30], the all-flavour neutrino spectrum from Cen A can be written as 

^ Cen A m = A- -m ( -M*) ^ Cen A (— ) > ( 2 - 14 ) 



sn'/i/n \'lvn J -'\>t '/i/n^br/ \>\vn / 

where £j (i = u, n) is the fraction of the proton's energy that is transferred to the species i 
in photopion interactions and t] vn is the ratio of the average neutrino energy to the average 
neutron energy. The KT model uses for these parameters the values featured in [30] , obtained 
from Monte Carlo simulations: f„ « 0.1, £ n « 0.5, (E u )/E p « 0.033 and (E n )/E p ps 0.83, 
with which £ v /£ n = 0.2 and rj un = 0.04. The neutrino break energy, E^, is estimated from 
the gamma-ray break energy as E^ r ~ 3 x 10 8 i? 7j br- Ref. [30] uses E 1 ^ 0Y = 200 MeV, so that 
i^br = 10 8 GeV. Under the assumption of equal flavour ratios at Earth, eq. (2.4), the + 
flux is 1/3 the flux in eq. (2.14). Plugging the power-law proton spectrum, eq. (2.5), with 
the normalisation constant for a point source, eq. (2.9), into the eq. (2.14) yields 

a (£ „, = *?- A Wfrf-'s^l (%.y (>L\ min L £.) (2 . 15) 

M 3 £ n £th \EthJ \-Ev,br/ V E V) t, T J 

for the muon-neutrino flux from Cen A, with E U ^ I = ri un E\> v = 4 x 10 6 GeV. Using the 
scaling relation, eq. (2.13), the muon-neutrino diffuse flux in the KT model is therefore 
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{E v ) ~ 5H (E th ) <j)^ cn A (E u ) and we can write it as 



C KT {E v ) = Af^T E - a min IK _K_ \ (2 16) 



with the neutrino normalisation given by 



"J sn 

and, following eq. (2.9), ^ cn A = $£ cn A ^ _ ^ 
2.3 Model by Becker & Biermann 

The BB model [16] describes the production of high-energy neutrinos in the relativistic jets 
of radio galaxies. According to the model, the UHECRs observed by the PAO originated at 
FR-I galaxies (relatively low-luminosity radio galaxies with extended radio jets, and radio 
knots distributed along them), which can in principle accelerate protons up to about 10 20 eV. 
Like in the KT model, here the protons are also shock- accelerated. Unlike the KT model, 
though, where the neutrino emission occurred in a region close to the AGN core, in the BB 
model the neutrino emission from p'j interactions is expected to peak at the first strong shock 
along the jet, lying at a distance of Zj ~ 3000 gravitational radii from the center. 

The optical depth corresponding to proton interactions with the disc photon field 
T pidisc ~ 0-02 an d so p'j interactions in the disc are not the dominant source of neutri- 
nos. The proton-proton interactions that occur when the jet encounters the AGN's torus are 
also neglected as neutrino source in the BB model. The dominant mechanism of neutrino 
production is the interaction between the accelerated protons and the synchrotron photons 
in the relativistic jet, at one of the jet's knots. For boost factors of the streaming plasma of 
T ~ 10, the optical depth T P7synch ~ 1. 

Hence, it is expected that neutrino emission occurs predominantly at the foot of the jet, 
where the beam is still highly collimated. Therefore, the BB model predicts a highly beamed 
neutrino emission, produced in the first shock {zj ~ 3000r s ), and consequently observable 
only from sources whose jets are directed towards Earth. Flat-spectrum radio sources, such as 
FR-I galaxies whose jets are pointing towards Earth, will have correlated neutrino and proton 
spectra, while steep-spectrum sources, which are AGN seen from the side, are expected to 
be weak neutrino sources, but to contribute to the cosmic-ray proton flux. 

The BB model assumes that the iV corr = 29 events that were observed by the PAO to 
have a positional correlation to sources in the VCV catalogue were indeed originated at AGN 
lying in the supergalactic plane. In order to relate the proton and neutrino normalisation 
constants, Ap lS and A„, we will use the connection between the proton and neutrino energy 
fluxes [16], i.e., 

_ r p7 T u 9, p n u 
3v ~ 12r CR ^n p { ^ CR)jp ' (2 ' 18) 

where Sl v , f^cR are the solid angles of emission of neutrinos and cosmic rays, respectively, 
and T u , Tqr are the boost factors of neutrinos and cosmic rays, respectively. The parameter 
z m |t x is the redshift of the farthest AGN that contribute to the cosmic-ray flux. The total 
number of neutrino (proton) sources, n„ (ra p ), is calculated by integrating the luminosity 
function of Willott [31] (Dunlop & Peacock [32]) from z 1 ^ = 0.018 (0.0008) up to zgg*. 
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On the other hand, assuming a power-law behaviour for the diffuse differential flux of 
protons, eq. (2.5), the energy flux results in 

jp = Ap f Ep — Ep ^E dEp = J A p (Op " 2)" 1 E vX 2 • if a P * 2 , (2.19) 

JE Pim in dE p y Ap In (-Ep, ma x/ -E'p.min) ; if %, = 2 

where the term proportional to -E p ,max 2 has been neglected, in the case when a p 7^ 2. As- 
suming that the neutrino spectrum follows the proton spectrum, i.e., <fit ' BB = A^E~ a " 
with ol u ~ atp, the energy flux for neutrinos is 

f («p - 2)" 1 , if a p + 2 

^ \^ iff ln(^ max /^,mm) , if Op = 2 l ' ; 

The lower integration limits for protons and neutrinos are, respectively, -E^min = F p m p ~ 
r p • (1 GeV) and E v>min = T u ■ (m T /4) = Y v ■ (0.035 GeV). Finally, replacing eq. (2.19), 
eq. (2.20), and the proton normalisation constant Ap s given by eq. (2.8) evaluated with 
-^evts = -^corD we see that when a.p 7^ 2, the neutrino normalisation constant is 

At m > BB - § (^)° P+1 ^ (^) QP ~ 2 < S (o„) • (2.21) 

The dependence of n u /n p on Zqr* is shown graphically in Figure 5 of Ref. [16]: n u /n p 
decreases with z^R*- To arrive at this expression 1 , it must be noted that because of the 
relativistic beaming in the jets, the emission solid angles are Q u ~ and Q p ~ 1/Tq R . 

When Q p = 2, the logarithms in the two spectra are similar and cancel out, making the 
previous expression for A^ lff ' BB valid also for a p = 2. Note that, since the ratio iV corr /H has 
decreased approximately by a factor of 2 between the original and updated PAO analyses, 
then the updated BB diffuse flux is about half the original. 



3 Current and preliminary bounds on the neutrino flux 

In the present work, we have assumed that the UHE AGN neutrino flux accounts for all 
of the UHE neutrino flux. This is, of course, a simplifying assumption, since high-energy 
contributions could also originate at other types of sources, such as gamma-ray bursts [33- 
36]. 

We have taken into account three experimental bounds on the diffuse astrophysical 
neutrino flux: two upper bounds, one set by the AMANDA-II experiment and the other by 
its successor, IceCube, in its half-completed configuration of 40 DOM strings; and a lower 
bound given by the discovery potential of the final 86-string IceCube configuration. These 
three bounds have been included in figure 1. 

The AMANDA-II upper bound on the diffuse high-energy flux of extra-terrestrial muon- 
neutrinos was obtained by using data recorded between the years 2000 and 2004 [38]: 

El<F^ < 7.4 x 10~ 8 GeV cm" 2 s" 1 sr^ 1 (90% C.L.) , (3.1) 

in the range 16 TeV - 2.5 PeV. This bound was set using exclusively upgoing UHE neutrinos, 
six of which were detected during the 807 days of live time reported. 

1 The reader should be wary that in their paper [16], Becker & Biermann incorrectly reported a dependence 
of the form ~ (r„/rcii) J-<lp due to an algebraic mistake[37]. 
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More recently, the IceCube Collaboration presented a preliminary upper bound using 
375 days of recorded upgoing data with the half-completed IceCube-40 array which is almost 
an order of magnitude tighter than the AMANDA bound [39]: 

El$f < 8 x 1(T 9 GeV cm" 2 s" 1 sr" 1 (90% C.L.) , (3.2) 

in the range 10 4 - 5 - 10 7 GeV. 

Finally, the discovery potential at the 5a level of the full, 86-string, IceCube array has 
been recently estimated [40] to reach, after five years of exposure, 

El(t>if < 7 x 10" 9 GeV cm" 2 s" 1 sr" 1 (5<r) , (3.3) 

also in the range 10 4 - 5 - 10 7 GeV. This is the estimated minimum necessary flux required 
for a 5<r discovery after five years of running IC86. We will use this discovery potential as 
a lower bound on the neutrino flux. The discovery potential in eq. (3.3) is better than the 
original estimate of 9.9 x 10 -9 GeV cm" 2 s" 1 sr" 1 that was presented in [41] due to a better 
knowledge of the detector and improved simulations. 

Note that these three bounds were obtained under the assumption of an E~ a neutrino 
flux, with q = 2. For the KT and BB models in our work, however, we have allowed for 
Therefore, we have calculated for each one of them the associated number of muon-neutrinos 
in the AMANDA, IceCube-40 and IceCube-86 configurations, as appropriate, by assuming 
an E~ 2 flux, and used these derived bounds on the number of events, and not on the flux, to 
constrain the KT and BB models. Concretely, we have assumed a cf>„ 1 ^ (E u ) = kE~ 2 flux, with 
the normalisation, k, given in each case by the numerical value of the bounds in eqs. (3.1)- 
-(3.3), in units of GeV -1 cm" 2 s" 1 sr . These numbers are displayed in table 1. The 
expressions required to calculate the number of upgoing muon-neutrinos in the AMANDA, 
IceCube-40, and IceCube-86 arrays, for an arbitrary diffuse neutrino flux 4>^, are contained 
in Appendix A. We have assumed that the effective detector area of AMANDA is 1 /100 times 
that of IceCube-86 and that the effective area of IceCube-40 is half the area of IceCube-86, 
on account of half the number of strings having been deployed. Note, however, that this is 
only an estimate, since the actual effective area of IceCube-40 will be strongly dependent on 
the efficiency of the cuts employed to calculate it. 

4 Muon-neutrino number of events in the IceCube-86 detector for the BB 
and KT models 

In this Section, we study the Koers &: Tinyakov (KT) and the Becker & Biermann (BB) mod- 
els of diffuse AGN neutrino flux through their predictions of the number of muon-neutrinos 
that will be detected by the full IceCube-86 neutrino detector. To calculate the number 
of neutrinos, we have adopted the method followed in [15], which is summarised here in 
Appendix A. 

In our analysis, we have fixed the IceCube-86 detector exposure time at T = 5 years and 
calculated the integrated event yield within the energy range 10 5 — 10 8 GeV. Only upgoing 
neutrinos have been considered, i.e., those that reach the detector with zenith angles between 
90° and 180° (the normal to the South Pole lies at 0°), for which the atmospheric neutrino 
and muon background is filtered out by interactions inside the Earth. Downgoing neutrinos, 
i.e., those with zenith angles between 0° and 90°, traverse only about 10 km of atmosphere 
before reaching the detector and have not been included in the analysis due to the added 
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Limit 


Energy range [GeV] 


Exp. time 


Upgoing 


AMANDA upper bound (AMANDA) [38] 
IccCube-40 preliminary upper bound (IC40) [39] 
IceCube-86 estimated 5cr discovery (IC86) [40] 


1.6 x 10 4 - 2.5 x 10 H 

io 4 - 5 - 10 7 
io 4 - 5 - 10 7 


807 days 
375 days 
5 years 


6.0 

5.90 
50.28 



Table 1. Maximum number of upgoing muon- neutrinos allowed by the reported exclusion limit from 
AMANDA and the preliminary one from IceCube-40, and minimum number of events needed for 
5cr discovery according to the estimated IccCube-86 5-year discovery flux. In every case, the event 
numbers were calculated by assuming a E~ 2 diffuse flux. Each bound on the number of events was 
calculated in the respective detector configuration (Appendix A contains the effective area for each), 
with the corresponding exposure time. 



Limit 


no source evolution 


strong source evolution 


AMANDA 


3.04 


2.81 


IC40 


2.59 


2.27 


IC86 


2.57 


2.25 



Table 2. Maximum value of the spectral index a in the Koers-Tinyakov model allowed by the upper 
bounds AMANDA and IC40, and minimum value needed for 5tr discovery according to the estimated 
IC86 discovery potential. 



difficulty of separating the atmospheric background from the astrophysical neutrino signal. 
Furthermore, in the case of the KT flux, we have considered both the scenario with no source 
evolution and the one with strong source evolution. 

Based on the experimental bounds introduced in the previous Section, we have defined 
two visibility criteria with the purpose of identifying the regions of parameter space allowed by 
the upper limits and accessible by the discovery potential of the full IceCube-86 array. Under 
the first one -the AMANDA visibility criterion-, the IceCube-86 event-rate predictions, for 
either KT or BB, are required to lie above the IC86 discovery potential and below the 
AMANDA upper bound. Similarly, under the IC40 visibility criterion, the event rates must 
lie above the IC86 discovery potential and below the IC40 upper bound. 

4.1 Parameters under study and neutrino fluxes 

We have calculated our expectations of the neutrino flux models taking as free parameters a 
for the KT model (to simplify, we will use a = a p hereafter), and a, T u /Tcr, and z™^* for 
the BB model, and varied them within the following intervals: 

2 < a < 3 , 1 < T u /T C r < 20 , 10 -3 < z^ < 0.03 . (4.1) 

This range of a has been chosen in order to cover a wide range around 2.7, the preferred value 
obtained from fits to combined cosmic-ray data [42], or values less than 2.3 that are predicted 
in case of stochastic shock acceleration [43-47]. We have defined the range of T^/Vcr for 
values greater than 1 since, under the assumptions made by the BB model, the neutrinos are 
produced in early shocks and protons, in late ones. Besides, it includes the value of 3 used in 
[16]. The range of z™R X ^ s ^ ne same as the one used in said reference. Our purpose in varying 
the latter parameter, z^R*' ^ s t° test different hypotheses about the maximum redshift up to 
which the AGN contribute to the UHE diffuse neutrino flux. We remind the reader that the 
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Figure 1. AGN muon- neutrino fluxes, multiplied by E 2 , according to the models by Becker & 
Biermann (BB) and Koers & Tinyakov (KT), with strong source evolution and without source evo- 
lution. The regions were generated by varying the model parameters in the ranges 2 < a < 3, 
1 < r„/r C R < 20, and 10~ 3 < zgg* < 0.03. The grey region corresponds to all the possible BB 
fluxes resulting from the variation of a, T v /Tcr, and z™R X ' whereas the brown and orange regions 
correspond to all the possible KT fluxes resulting from the variation of a, under the assumption of 
no source evolution and of strong source evolution, respectively. The atmospheric muon-neutrino 
flux has been plotted (in black, dotted, lines) for comparison. The AMANDA-II upper bound, the 
preliminary 40-string IceCube upper bound and an estimated 86-string IceCube five-year discovery 
potential at 5cr have been included by assuming a E~ 2 flux. The atmospheric neutrino flux is given 
by the parametrisation in Ref. [15]. See the text for details. 



results for the KT model have been obtained for a fixed value of -Zqr x = 5 and so they were 
not affected by this variation. 

Figure 1 shows the BB and KT diffuse muon-neutrino fluxes, multiplied by E^,, as 
functions of the neutrino energy, when the values of the model parameters are varied within 
the ranges that we have quoted above. We have also included the upper bounds on the 
flux set by AMANDA and IceCube-40, and the estimated discovery potential of IceCube-86 
after five years of running. Our analysis will focus on the different regions enclosed between 
these upper bounds and the IC86 discovery potential taken as a lower bound, in the energy 
range 10 5 -10 8 GeV, where the fluxes may be detected in IceCube. We will find how the 
bounds on the neutrino flux translate into bounds on the values of a, T u /Tcr, and ^™R X ' 
thus restricting the capacity of the KT and BB flux models to account for an observed 
extra-terrestrial neutrino signal. 

4.2 KT event-rate expectations in IceCube-86 

Since the KT flux depends on a single parameter, i.e., the spectral index a, we can translate 
the bounds on event numbers directly into bounds on a. In this way, the results presented 
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Figure 2. Integrated number of upgoing muon- neutrinos, between 10 5 and 10 s GeV, expected in 
IceCubc-86, after T = 5 years of exposure, associated to the KT production model assuming (a) no 
source evolution and (b) strong source evolution. The orange-coloured bands are the regions of values 
of a that lies above the IC86 discovery potential and below the IC40 upper bound, while the hatched 
region lies above IC86 and below the AMANDA upper bound (see table 2). 



in table 2 represent the upper limits on a given by the AMANDA and IC40 bounds and the 
lower limits given by the IC86 discovery potential. 

Figure 2 shows the integrated number of upgoing muon-neutrinos with energies between 
10 5 and 10 8 GeV, as a function of a, that is expected in the full 86-string IceCube array after 
five years of exposure. Plot (a) assumes no source evolution, whereas (b) assumes strong 
source evolution. The predictions under the assumption of strong source evolution are up 
to an order of magnitude higher than under no source evolution. This fact can be easily 
understood since a difference of a similar magnitude is found in the neutrino boost factor, as 
shown in Ref. [15]. 

The orange-coloured and hatched bands mark the visibility regions under the IC40 and 
AMANDA visibility criteria, respectively, according to table 2. Owing to the fact that the 
AMANDA upper bound is less restrictive than the IC40 bound, the visibility regions are in 
every case larger when the former one is used. According to figure 2 and table 2, the ranges 
of event numbers, A^kis that IceCube-86 will be able to detect in the interval 10 5 — 10 s GeV, 
after five years of exposure, are: 

68 < A£t < 77(1847) , (4.2) 
assuming no source evolution and using the IC40 (AMANDA) upper bound, and 

85 < < 95 (2709) , (4.3) 

assuming strong source evolution. 

From table 2, we see that the KT model with no source evolution is allowed for higher 
values of a than the model with strong source evolution. This is due to the fact that the 
KT flux grows with a, and that, for a given value of a, the event yield produced by the 
strong source evolution model is up to an order of magnitude higher than the yield with no 
source evolution. Thus, lower values of a are needed to keep the former below the IC40 or 
AMANDA event-number upper bounds. 

From the same table, we find that for the KT model the value of a = 2.7, obtained from 
fits to cosmic-ray data, would still be allowed under the AMANDA visibility criterion, but 
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Figure 3. Variation of the integrated number of upgoing muon- neutrinos expected in the range 
10 5 < Ev/GcV < 10 8 associated to the BB model, after T = 5 years of exposure of the IceCube-86 
detector. In (a), (b), and (c), the value of z^r* nas been fixed, respectively, at the representative 
values of 10~ 3 , 0.01, and 0.03, while a and T v /Tcr have been allowed to vary. Likewise, in (d), (e), 
and (f), r^/TcR has been fixed at 1, 10, and 20, respectively, while a and z™^* have been varied. 
The solid lines are iso-contours of number of events: 10 (solid black), 100 (dashed red), 300 (dotted 
blue), and 500 (dash-dotted green). The region coloured orange is the parameter region where the 
event-number predictions lie above the IC86 discovery potential and below the IC40 upper bound, 
i.e., the IC40 visibility region. Similarly, the hatched region is where the predictions lie above the 
IC86 potential and below the AMANDA upper bound, i.e., the AMANDA visibility region. 



is discarded by the more recent IC40 criterion, regardless of the choice of source evolution. 
Under the assumption of strong source evolution, the other proposed value of a = 2.3 is 
excluded (permitted) by the IC40 (AMANDA) visibility criterion, while values of a < 2.25 
would be out of reach of the IceCube discovery potential. Under no source evolution, the 
region of a below the IC86 potential starts from 2.57. This constitutes a strong hint toward 
the KT flux being too large. However, as explained in Section 3, we would like to stress 
that our visibility criteria make use of event-yield bounds that are deduced from bounds on 
a E~ 2 flux, a comparison that might be overly reducing the size of the visibility regions. A 
more sophisticated analysis that makes use of model- independent flux bounds, i.e., bounds 
not exclusive to E~ 2 models, will be presented elsewhere [48]. 

4.3 BB event-rate expectations in IceCube-86 

As to the BB flux model, figure 3 shows iso-contours of the expected integrated number of 
upgoing muon-neutrinos in the IceCube-86 detector, in the T v /Tqr-oi plane, for fixed values 
of (a) zgp* = 10~ 3 , (b) 0.01, and (c) 0.03, and in the z^^-a plane, for fixed values of (d) 
r^/rcR = 1, (e) 10, and (f) 20. The BB normalisation constant, according to eq. (2.21), 
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Z CR 


Minimum 


Maximum a 


Minimum 


Maximum r, y /rcR 




a 


IC40 visib. 


AMANDA visib. 


r,/r CR 


IC40 visib. 


AMANDA visib. 


10" a 


2 


2.3 


2.65 


l 


3 


7.5 


0.01 


2 


2.6 


2.95 


l 


8.5 


20 


0.03 


2 


2.65 


3 


l 


11 


20 



Table 3. Allowed intervals of a and T^/Tcb. obtained by projecting the visibility regions from plots 
3a-c onto the axes. 



r^/rcR 


Minimum 

a 


Maximum a 


Minimum zgp x 


Maximum 

-max 
Z CR 


IC40 visib. 


AMANDA visib. 


IC40 visib. 


AMANDA visib. 


i 


2.25 


2.65 


2.9 


10" a 


lO" 3 


0.03 


10 


2 


2.03 


2.22 


0.015 


0.002 


0.03 


20 


2 


2.03 


2.1 


0.015 


0.008 


0.03 



Table 4. Allowed intervals of a and zgp x obtained by projecting the visibility regions from plots 3d-f 
onto the axes. 



decreases with Zqr* and increases with T u /Tqr- This behaviour is observed in figure 3, where, 
for fixed values of a and T u /Tcr, the number of events decreases as ^§R X increases. On the 
other hand, for fixed values of a and z§r x , the number of events increases with T u /Tcr- 

In each plot, as we have mentioned before, the IC40 visibility region is coloured orange 
and lies between the IC86 discovery potential (left border) and the IC40 upper bound (right 
border) listed in table 1. The AMANDA visibility region, on the other hand, is represented 
by the hatched region, and its right border is fixed instead by the AMANDA bound. 

Besides the observed narrowness of the visibility regions, there are two main features 
to point out. First, if the value of Zq^ increases, the allowed ranges of a and T u /Tcr also 
increase, with higher values being allowed. Second, if the value of T u /Tqr increases, the 
allowed ranges of a and z^R* decrease, with a tending to lower values and 2™R X to higher 
ones. These observations can be quantified if we project the visibility regions in each plane 
onto the horizontal and vertical axes. The allowed regions of the parameters are shown in 
tables 3 and 4. 

In light of the results presented in these tables, and momentarily assuming that a = 2.7 
is the true value of the cosmic-ray spectral index [42], we see that under the AMANDA 
visibility criterion the BB flux model is clearly excluded for T u /Tqr > 10 (for any value 
of Zq^) and also for the lowest values of £™]| x , close to 10~ 3 (for any value of T u /Tqr)- 
Whenever a = 2.7 is allowed by the BB model, it is only inside a very narrow region of 
parameter space, around IV/Tcr ~ 1 and ^§R X ~ 0.004. On the other hand, under the more 
recent IC40 visibility criterion, the BB model at a = 2.7 is discarded for all values of T u /Tqr 
and zgg*. 

If we consider the other values of a = 2.3 and 2.0 proposed in the literature (see Ref. [16] 
and references therein), we find that the allowed regions, for a = 2.3 and the AMANDA 
visibility criterion, are: 1 < IY/r C R < 3, 2.5 < T v /T C r < 6.5 and 3 < IV/Tcr < 8 for 
z gg* = 10" 3 ,0.01 and 0.03, respectively. In the case of a = 2, and the AMANDA visibility 
criterion, the allowed regions are: 3 < T u /T C r < 8, 8 < T^/r C R < 20 and 11 < T V /T CR < 20 
for zg^ x = 10~ 3 ,0.01 and 0.03, respectively. For a = 2.3(2.0), and the IC40 visibility 
criterion, the allowed values for T u /T C r are: 1(3), 2.5(8) and 3(11) for z^ x = 10~ 3 , 0.01 and 



- 13 - 




Figure 4. Separation between the BB and KT models, in terms of A = |7Vbb — Nkt\, measured in 
units of a = \/ Nkt (see text), for upgoing neutrinos with energies in the range 10 5 < E u /GeV < 10 8 
and assuming no source evolution for the KT model. The exposure time T = 5 years. In (a), (b), and 
(c), the value of .z™]| x has been fixed, respectively, at the representative values of 10~ 3 , 0.01, and 0.03, 
while a and T^/Tck have been allowed to vary. Likewise, in (d), (e), and (f), T u /Tqr has been fixed 
at 1, 10, and 20, respectively, while a and z™^* have been varied. The solid lines are iso-contours of 
A = la (solid black), 5a (dashed red), 10cr (dotted blue), and 20<r (dash-dotted green). The region 
coloured orange is the parameter region where the event-number predictions of, simultaneously, the 
KT and BB models lie above the IC86 discovery potential and below the IC40 upper bound, i.e., 
the IC40 visibility region. Similarly, the hatched region is where the predictions of both models lie 
above the IC86 potential and below the AMANDA upper bound, i.e., the AMANDA visibility region. 
These regions of simultaneous visibility are where comparison between the two production models is 
meaningful, according to each of the two visibility criteria. 



0.03, respectively. Clearly, lower values of a fare better under the more recent IC40 upper 
bound. Like for the KT model, the BB model region of parameter space could be larger if 
an analysis based on non-E~ 2 bounds were performed instead. 

5 Comparison between the KT and BB models using the IceCube detector 

We have quantified the difference between the predictions put forward by the two models 
using the quantity 

A(a,iyr C R,sgr) = I^VbbK Wr CR ,z^ x ) -iVKxHl , (5.1) 

and expressed it in units of a (a) = \J Nkt (a), he., at every point in parameter space we 
have measured the difference between the number of events predicted by each model, in units 
of the standard deviation of the KT prediction, assuming for it an uncertainty characteristic 



- 14 - 




Figure 5. Same as figure 4, but assuming strong source evolution for the KT model. 



of a Gaussian distribution. The higher the value of A, the greater the difference between 
the predictions. The comparison between the models, however, is only valid within the 
region that results from the intersection of the individual KT and BB visibility regions, 
given, respectively, by table 2 and figure 3. This guarantees that the numbers of events 
predicted by both models lie above the minimum required signal for detection at 5a from the 
atmospheric neutrino background, so that the comparison between them is meaningful. 

Figures 4 and 5 show the separation between the models using the integrated number of 
muon-neutrinos in the IceCube-86 detector. The iso-contours correspond to A ja = 1 (solid 
black), 5 (dashed red), 10 (dotted blue), and 20 (dash-dotted green), in the plane r„/TcR- 
-a, for values of (a) zg^ x = 10~ 3 , (b) 0.01, and (c) 0.03, and in the plane log (z^)-a, 
for values of (d) T u /Tcr = 1, (e) 10, and (f) 20. Where only one or none of the models 
are visible, the discrimination between them is obvious or meaningless, respectively. We 
have coloured orange the region of simultaneous visibility under the IC40 criterion, and 
hatched the region of simultaneous visibility under the AMANDA criterion. Evidently, since 
the individual visibility regions of the KT and BB models are larger under the AMANDA 
visibility criterion than under the IC40 criterion, the regions of simultaneous visibility are in 
every case larger under the former. 

We see that the KT and BB visibility regions overlap only at low values of T u /Tqh 
and that the size of the overlapping regions grows with z™^ x , so that they are largest for 
Z CR X = ^-01 and 0.03, as shown in plots (b) and (c) of figures 4 and 5. In particular, 
under the AMANDA visibility criterion, and assuming no source evolution, the regions of 
simultaneous visibility exist only for low values of T u /Tcr, between 1 and 3, while assuming 
strong source evolution, they exist up to r„/TcR ~ 10. Under the IC40 visibility criterion, 



- 15 - 



comparison is allowed only inside very small regions of simultaneous visibility that lie at 
a ~ 2.57(2.25) - 2.59(2.27), r„/r CR ~ 1(3) - 1.5(4), and zgg* = 0.01 - 0.03, assuming 
no (strong) source evolution. Hence, comparison between the models becomes unfeasible in 
most of the parameter space. 

Regardless, within the small IC40 simultaneous visibility region, the models can be 
separated in no less than 5a and no more than 10c, under both assumptions on source evo- 
lution, whereas under the dated AMANDA visibility criterion separations can vary between 
la and 20<r. Separations of 5a would be sufficient to discern in a statistically meaningful 
way between the KT and BB models. Notice that the comparison at the favoured value of 
a = 2.7 is not allowed under the IC40 visibility criterion, since neither flux will be visible 
in IceCube-86. For a = 2.0 and 2.3, there is no region of simultaneous visibility under this 
same visibility criterion. 

6 Summary and conclusions 

We have studied the IceCube-86 event rate expectations for two models of AGN diffuse 
muon-neutrino flux proposed in the literature, one by Koers & Tinyakov (KT) [15] and 
another by Becker & Biermann (BB) [16], both of which take into account the apparent 
correlation, reported by the Pierre Auger Collaboration [7], between the incoming directions 
of the highest-energy (E > 55 EeV) cosmic rays and the positions of AGN in the 12th 
edition Veron-Cetty & Veron catalogue [9]. In doing this, we have assumed that the flux 
of neutrinos from AGN makes up all of the UHE astrophysical neutrino flux. Both models 
propose a power-law flux, i.e., proportional to E~ a , resulting from shock acceleration. 

In our analysis, we have taken the spectral index, a, as well as two other parameters 
associated to the BB model, namely, the ratio of relativistic boost factors of neutrinos and 
cosmic rays, T u /Tqr, and the redshift of the most distant AGN that contributes to the 
diffuse cosmic-ray flux, z™r x , as free parameters, and varied their values within the following 
intervals: 2 < a < 3, 1 < T u /T C r < 20, and 10~ 3 < zgg* < 0.03. In addition, we have 
explored the KT model under two assumptions on the evolution of the number density of 
AGN: either they do not evolve with redshift, or they evolve strongly with it, following the 
star formation rate. Neutrino fluxes calculated using the latter assumption are up to an order 
of magnitude higher than the ones calculated using the former one. 

For each point (a, T u /Tcr, Zq^) in parameter space, we have calculated for both models 
the associated integrated number of upgoing muon-neutrinos, between 10 5 and 10 8 GeV, 
that is expected after five years of exposure of the full 86-string IceCube neutrino detector 
(IceCube-86). In order to determine the regions of parameter space that this detector will 
be able to probe, we have tested two different upper bounds on the UHE neutrino flux: the 
bound reported by the AMANDA Collaboration using 807 days of observation [38] and a 
preliminary bound obtained after 375 days of exposure of the half-completed IceCube-40 
detector (IC40) [39]. A lower bound, on the other hand, was fixed at the estimated IceCube- 
86 five-year discovery potential at the 5a level (IC86) [40]. With this we have defined "regions 
of visibility" in parameter space as those regions inside which the event-rate predictions lie 
above the IC86 discovery potential and below the AMANDA or IC40 upper bound. Since 
the IC40 upper bound is lower than the AMANDA bound, the former restricts the allowed 
parameter space more than the latter. 

It is possible to confine the spectral index of the KT model within the range 2.57 < a < 
2.59(3.04), under the assumption of no source evolution and using the IC40 (AMANDA) 
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upper bound, and 2.25 < a < 2.27(2.81), under the assumption of strong source evolution. 
For the BB model, we found that IceCube-86 is sensitive to high values of T u /Tcr, close 
to 20, only within small regions of parameter space, with a < 2.1 and £§R X ~ 0-03. For 
1 < T v /Tqr < 11, under the IC40 visibility criterion, the spectral index can take on values 
within the interval 2 < a < 2.65, though the highest values are accessible only with -z™^ x = 
0.01 to 0.03. For low values of T u /Tcr, around 1, the allowed ranges are 2.25 < a < 2.65 
and 10~ 3 < zg^ x < 0.03. 

Using combined cosmic-ray data [42], the preferred value of a has been set at 2.7. We 
have found that, if the AMANDA upper bound is used, this value is allowed in both the KT 
and BB models, whereas if the more recent IC40 upper bound is used, it is not. The authors 
of [16] claim that the true value of the spectral index might be either a = 2.0 or 2.3. For the 
BB model, these two values are allowed under both visibility criteria. For the KT model, 
using the AMANDA bound, the value a = 2.3 is allowed under strong source evolution, while 
under no source evolution it is not testable since it lies below the IC86 discovery potential. 
Using the IC40 bound, a = 2.3 is excluded under strong source evolution and is also not 
testable under no source evolution. The value a = 2.0 is not testable under any assumption 
on the source evolution. Note, however, that the experimental discovery potential and upper 
bounds that we have used were calculated for a E~ 2 flux and that using them to constrain 
the BB and KT models might be slightly over-constraining the parameter space. 

Additionally, in the event that an UHE neutrino signal is detected after five years 
of running the full IceCube array, and assuming that it was produced solely by the neu- 
trino flux from AGN, we have explored the detector's capability to distinguish between 
the KT model, with strong and no source evolution, and the BB model, i.e., to determine 
which one of the two models would correctly describe the detected UHE neutrino data. 
In order to do this, we have defined a measure of the separation between the models as 
A(a,T u /T C R,z^ x ) = \N BB (a, T u /T C r, z^ x ) - N KT (a)\, with A BB and A KT the number 
of muon-neutrinos expected in IceCube-86 associated to each model, between 10 5 and 10 8 
GeV, after five years of running. At each point in parameter space, we have calculated the 
value of A, expressed in units of a (a) = \J Akt (a). The comparison between the flux 
models, however, is meaningful only in those regions of parameter space where both models 
simultaneously lie inside their respective visibility regions. Thus, under the IC40 visibility 
criterion, comparison is allowed only inside very small regions of simultaneous visibility lo- 
cated at a ~ 2.57(2.25)-2.59(2.27), T u /T C r 1(3)-1.5(4), and = 0.01-0.03 assuming 
no (strong) source evolution. Within these regions, the separation between models is at the 
level of 5a or higher. Hence, comparison between the models becomes unfeasible in most of 
the parameter space, but where it becomes possible, it is statistically meaningful. 

A comment is in order: if, for the BB model, we had performed the integration in 
Z CR X U P t° a value < 5, the associated number of events would have been larger and the 
corresponding visibility region even tighter than the ones we have presented, for which the 
contributions to the diffuse flux only come from the supergalactic plane (z^^ < 0.03). 
Since the magnitude of the separation between models relies on the number of events, then 
either the level of separation would have been higher or there would have been no region of 
simultaneous visibility. 

We have thus shown that, after five years of running, the completed IceCube array 
might be able to strongly constrain the KT and BB models, leaving only small regions of 
parameter space where the models survive. In addition, discrimination between the models, 
while feasible only within even smaller regions of parameter space, might be able to reach 
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the 5<7 level. The reader should be aware that our predictions are based on an all-proton 
cosmic-ray flux, but there is growing evidence that the UHECR flux is composed mainly of 
heavy nuclei [49-52] (see, however, [53, 54]), and, as a consequence, the UHE neutrino flux 
would be reduced. Thus, with reservations, our results might be seen as symptoms of the 
need for new models of AGN neutrino production that are better equipped to face the latest 
experimental bounds on the UHE neutrino flux. 
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A Neutrino detection in IceCube 

We have calculated the predicted number of muon-neutrinos detected in IceCube-86 using the 
method presented in Ref. [15]. In general, the integrated number of upgoing muon-neutrinos 
at a Cerenkov detector due to a diffuse flux of muon-neutrinos, (/>^ ff , with energies between 
E™ m and £™ ax , is calculated as 

gmax 

N„ tVp = TSl " dE v <pf^ (E„) A* p eS (E v ) , (A.l) 

JE™ in 

where T is the detector's exposure time; f2, the detector's opening solid angle; E u , the 
neutrino energy; 0" 1 is either the KT or BB diffuse AGN neutrino flux; and A^ p cfS: is the 
upgoing neutrino effective area. 

Note that the six extra DeepCore strings of the IceCube-86 array increase the neutrino 
effective area only in the range 10 < E u /GeV < 10 3 [55]. Above 10 3 GeV, the IceCube 
effective area is determined solely by the remaining 80 strings. 

The effective neutrino area takes the form 

A Zs = S (E v ) P, A^ cS (EJ , (A.2) 

where S is the shadowing factor, which takes into account neutrino interactions within the 
Earth; P^, the probability that the neutrino-spawned muon reaches the detector with energy 
greater than the threshold energy E™ m required to be detected; and A^g, the detector's 
effective area for muons. We will explain each term in eq. (A.2) in what follows. 
The probability of muon detection can be written as [15] 

P M (£7 M ) = 1 - exp(-N Av a^ (E u ) , (A.3) 
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Figure 6. Angle-averaged upgoing neutrino effective areas, as functions of the neutrino energy, for 
the IceCube-86, IceCube-40, and AMANDA detectors. The IceCube-40 effective area is estimated at 
half the IceCube-86 area (see text), while the AMANDA effective area is a factor of 100 lower than 
the IceCube-86 area. 



where Aa v = 6.022 x 10 23 mol" 1 = 6.022 x 10 23 cm" 3 (w.e., water equivalent) is Avogadro's 
constant; is the charged-current neutrino- nucleon cross section, taken from [56] (which 
uses CTEQ4 data); and i? M is the muon range within which the muon energy reaches the 
threshold energy E™ m = 100 GeV, which can be expressed as 



M ^ b \a + bE™ m 



(A.4) 



with a = 2.0 x 10 -3 GeV cm -1 (w.e.) accounting for ionisation losses and b = 3.9 x 10~ 6 
cm" 1 (w.e.) accounting for radiation losses. The relation between neutrino and muon energy 
is obtained by assuming single-muon production in each neutrino interaction, which leads to 
E[i = ycc (E U )E U , with ycc the mean charged-current inelasticity parameter tabulated in 
[56]. 

The shadowing factor, S, is defined in terms of P v (E v , 9), the probability that a neutrino 
arriving at Earth with nadir angle 9 (the North Pole is located at 9 = 0°) and interacting 
with Earth matter, reaches the detector. We use [15] 



S(E U ) 



1 - COS (0 max ) Jo 



d9 sin (9)P V (E U , 9) , 



where (9 max is the detector's maximum viewing angle, which we have taken to be 9 V 
as in Ref. [15]. Thus, the detector's opening angle is 



(A.5) 

= 85°, 



n 



2- 



sin (6») d9 = 2vr [1 - cos (0 max )] « 5.736 sr 
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(A.6) 



where a^°jy is the total (charged- plus neutral-current) neutrino-nucleon cross section, tab- 
ulated in Ref. [56]; p(r) is the Earth's density profile given by the Preliminary Reference 



te = 6371 km the Earth radius; and L(6) = 2r£?cos(#) is the distance that a neutrino 
traversing the Earth at angle 9 propagates. 

Lastly, for IceCube-86's upgoing muon effective area, A^q, we have used the curve 
corresponding to level-2 cuts in Figure 5 of Ref. [41], which is the effective area averaged over 
the northern hemisphere, and dependent only on the incoming muon energy, E^. Figure 6 
shows that the IceCube-40 neutrino effective area is estimated at one half the IceCube-86 
effective area (see Section 3), while the AMANDA neutrino effective area was a factor of 100 
lower than IceCube-86 area. 
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